home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2002 November / SGI IRIX Base Documentation 2002 November.iso / usr / share / catman / p_man / cat3 / SCSL / ztpcon.z / ztpcon
Encoding:
Text File  |  2002-10-03  |  4.0 KB  |  133 lines
  1.  
  2.  
  3.  
  4. ZZZZTTTTPPPPCCCCOOOONNNN((((3333SSSS))))                                                          ZZZZTTTTPPPPCCCCOOOONNNN((((3333SSSS))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      ZTPCON - estimate the reciprocal of the condition number of a packed
  10.      triangular matrix A, in either the 1-norm or the infinity-norm
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE ZTPCON( NORM, UPLO, DIAG, N, AP, RCOND, WORK, RWORK, INFO )
  14.  
  15.          CHARACTER      DIAG, NORM, UPLO
  16.  
  17.          INTEGER        INFO, N
  18.  
  19.          DOUBLE         PRECISION RCOND
  20.  
  21.          DOUBLE         PRECISION RWORK( * )
  22.  
  23.          COMPLEX*16     AP( * ), WORK( * )
  24.  
  25. IIIIMMMMPPPPLLLLEEEEMMMMEEEENNNNTTTTAAAATTTTIIIIOOOONNNN
  26.      These routines are part of the SCSL Scientific Library and can be loaded
  27.      using either the -lscs or the -lscs_mp option.  The -lscs_mp option
  28.      directs the linker to use the multi-processor version of the library.
  29.  
  30.      When linking to SCSL with -lscs or -lscs_mp, the default integer size is
  31.      4 bytes (32 bits). Another version of SCSL is available in which integers
  32.      are 8 bytes (64 bits).  This version allows the user access to larger
  33.      memory sizes and helps when porting legacy Cray codes.  It can be loaded
  34.      by using the -lscs_i8 option or the -lscs_i8_mp option. A program may use
  35.      only one of the two versions; 4-byte integer and 8-byte integer library
  36.      calls cannot be mixed.
  37.  
  38. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  39.      ZTPCON estimates the reciprocal of the condition number of a packed
  40.      triangular matrix A, in either the 1-norm or the infinity-norm. The norm
  41.      of A is computed and an estimate is obtained for norm(inv(A)), then the
  42.      reciprocal of the condition number is computed as
  43.         RCOND = 1 / ( norm(A) * norm(inv(A)) ).
  44.  
  45.  
  46. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  47.      NORM    (input) CHARACTER*1
  48.              Specifies whether the 1-norm condition number or the infinity-
  49.              norm condition number is required:
  50.              = '1' or 'O':  1-norm;
  51.              = 'I':         Infinity-norm.
  52.  
  53.      UPLO    (input) CHARACTER*1
  54.              = 'U':  A is upper triangular;
  55.              = 'L':  A is lower triangular.
  56.  
  57.      DIAG    (input) CHARACTER*1
  58.              = 'N':  A is non-unit triangular;
  59.              = 'U':  A is unit triangular.
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. ZZZZTTTTPPPPCCCCOOOONNNN((((3333SSSS))))                                                          ZZZZTTTTPPPPCCCCOOOONNNN((((3333SSSS))))
  71.  
  72.  
  73.  
  74.      N       (input) INTEGER
  75.              The order of the matrix A.  N >= 0.
  76.  
  77.      AP      (input) COMPLEX*16 array, dimension (N*(N+1)/2)
  78.              The upper or lower triangular matrix A, packed columnwise in a
  79.              linear array.  The j-th column of A is stored in the array AP as
  80.              follows:  if UPLO = 'U', AP(i + (j-1)*j/2) = A(i,j) for 1<=i<=j;
  81.              if UPLO = 'L', AP(i + (j-1)*(2n-j)/2) = A(i,j) for j<=i<=n.  If
  82.              DIAG = 'U', the diagonal elements of A are not referenced and are
  83.              assumed to be 1.
  84.  
  85.      RCOND   (output) DOUBLE PRECISION
  86.              The reciprocal of the condition number of the matrix A, computed
  87.              as RCOND = 1/(norm(A) * norm(inv(A))).
  88.  
  89.      WORK    (workspace) COMPLEX*16 array, dimension (2*N)
  90.  
  91.      RWORK   (workspace) DOUBLE PRECISION array, dimension (N)
  92.  
  93.      INFO    (output) INTEGER
  94.              = 0:  successful exit
  95.              < 0:  if INFO = -i, the i-th argument had an illegal value
  96.  
  97. SSSSEEEEEEEE AAAALLLLSSSSOOOO
  98.      INTRO_LAPACK(3S), INTRO_SCSL(3S)
  99.  
  100.      This man page is available only online.
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.